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Abstract This paper aims to solve a class of CEC benchmark constrained 
optimization problems that have been widely studied by nature-inspired opti- 
mization algorithms. Global optimality condition based on canonical duality 
theory is derived. Integrating the dual solutions with the KKT conditions, we 
are able to obtain the approximate solutions or global solutions easily. 
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1 Introduction 

Nature-inspired optimization algorithms, such as genetic algorithm (GA), evo- 
lution strategy (ES), particle swarm optimization (PSO) and differential evo- 
lution (DE), have received considerable attention in recent decades due to 
their strong adaptability and easy implementation. Strictly speaking, these 
algorithms are unconstrained optimization procedures, and therefore it is nec- 
essary to find techniques to deal with the constraints when solving constrained 
optimization problems. The most common approach to handle constraints is 
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the penalty function method. The idea of this approach is to transform a con- 
strained optimization problem into an unconstrained one by adding a certain 
term to the objective function based on the amount of constraint violation. 
Then, some special representations and operators are designed to preserve the 
feasibility of solutions at all times or to repair a solution when it is infeasible. 
Multiobjective optimization techniques are also applied to manage constraints. 
The main idea is to rewrite the single objective optimization problem as a mul- 
tiobjective optimization problem in which the constraints in original problem 
are treated as additional objectives [51ITU]. 

By introducing a Lagrangian multiplier vector to relax the constraints, the 
classical Lagrangian is a saddle function if the objective function and all of the 
constraints are convex. Under certain constraint qualifications, the well-known 
strong min-max duality relation holds, and in this case, the problem can be 
easily solved by well-developed convex programming techniques. However, due 
to nonconvexity of either objective function or the constraints, the Lagrangian 
is no longer a saddle function and only the weak duality relation holds, leading 
to the duality gap in global optimization [I][2JE]. In order to bridge the gap 
inherent in the classical Lagrange duality theory, the canonical duality theory 
has been developed recently. Its core is to transform a nonconvex primal min- 
imization problem to the concave canonical dual maximization problem over 
a convex space without duality gap by a canonical dual transformation [7]. 

As shown by the global optimality condition contained in the canonical 
duality theory, if the dual solution is in the positive definite domain, it is easy 
to get the corresponding global solution to the primal problem. However, if the 
condition is not satisfied, some strategies are necessary to recover the global 
solution [51I51IT2"]. In this paper, we focus on solving a class of CEC (Congress 
on Evolutionary Computation) benchmark constrained optimization problems 
that have been widely studied by nature-inspired algorithms. By integrating 
the canonical duality theory with the KKT conditions, we are able to obtain 
the approximate solutions or global solutions easily. 



2 The canonical duality theory 

In this paper, we focus on the following quadratic optimization problem with 
quadratic and box constraints (primal problem): 

(V) : min |p(x) = ^x r Ax - a T x - a : x 6 K n j, 

s.t. g(x) = {.9j(x)} = {^Bjx-bjx- bj} < 0,j = 1, ••• ,m, 
d < Xi < di,i = 1, • • • , n, (1) 

where, x = (x±, ■ ■ ■ ,x n ),A — A T , Bj = Bj g IR nx ™ are symmetric matrices, 
a, bj e R™ are given vectors, a, bj, Ci, di are constant. 

Let Ei £ M. nxn , Bi £ K™ be a diagonal matrix and a unit vector, with 
all zeros except a one in the position (i,i) and (i), respectively, Bk = 2Ek, 
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bfe = (cfe + dk)ek, bk = Ckdk, k = m + 1, • • ■ ,m + n, then constraints in (V) 
can be uniformly rewritten to 

g(x) = { 3fe (x)} = {ix T B fc x - b£x - b k } < 0, k = 1, • • • ,m + n. (2) 
Firstly, we introduce an indicator function: 

f if e < 

W(e) = \ - . (3) 

I +oo otherwise 

where, e = (ei, • • • , e m +„), and let C/(x) = — /(x) = — ix T Ax + a T x + a, then 
the primal problem [V) can be written in the following form: 

min{n (x) = W(g(x)) - C/(x) : x G W 1 }. (4) 

Secondly, we introduce a nonlinear operator 

i = (tu ■ ■ ■ ,U+n) = {6} - A(x) = {ix T B fc x - b£x - b k } (5) 

so that g(x) can be recast by V(£,) — £ = g(x), and then the primal problem 
(V) can be reformed in the canonical form: 

min{n(x) = W(V(A(x))) - U(x) : x e R"}. (6) 

By the Fenchel transformation, the conjugate function W^(cr) of W(e) can 
be defined by 

„ rp { if er > 

W^a)= S up{e T a-W(e)}={ (7) 
e I +oo otherwise 

which is convex and lower semi-continuous. 

According to the relations er e dW(e) «e£ dW*(cr) ^ W(e) + VF J (er) = 
e T cr from convex analysis, we can replace VF(g(x)) by g T (x)er — VF B (er), and 
then we can get the extended Lagrangian 

~o(x, tr) = g T (x)<r - wHa) - [/(x). (8) 

Next, we introduce the invertible duality mapping 

C=(?l,"-, Wn) = {ft} = W(£) = /• (9) 

Defining the Legendre conjugate 1^*(<0 = sta{£ T £ — V(£,)}, and using the 
equivalent relations ^ = W(£) £ = VV*(s) 4^ £ T <j = V(£)+V*(s), we can 
replace g(x) by A T (x)<; — V*(s), so we obtain the generalized complementary 
function 

~(x, cr) = ix T G(cr)x - x T F(er) - cr T d - a, (10) 
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where. 

G(o-) = A + <7kB k ,F((r) = a + ^ °kb k , 

k=l k=l 

and 

d = (&!,••■ ,U)V-K-,^+»)eiT + " = {^K m+ >>0}. 

By using the generalized complementary function, the canonical dual func- 
tion P d (cr) can be formulated as 

P d (cr) = sta{~(x,<7)}. (11) 

Solving the critical points of H(x, cr), we can get the canonical equilibrium 
equation 

G(<r)x = F{&). (12) 

For any given cr, if F{cr) is in the column space of G(er), denoted by 
C j(G(<x)), i.e., a linear space spanned by the columns of G(er), the solution 
of the canonical equilibrium equation can be well defined by 

x = G t (<x)F(<x), (13) 

where, G'(<j) denotes the Moore-Penrose generalized inverse of G(cr). 
Then, the canonical dual function can be written explicitly as follows 

P d (cr) = -^F t (ct)G\<t)F((t) - a T d - a. (14) 
Finally, the canonical dual problem can be posed as follows: 

(V d ) : max{PV) ■ <r € 5+} (15) 
where, the dual feasible space is defined by S+ = {cr G R™ + "|G(er) y 0}. 

Theorem 1 (Global Optimality Condition). Suppose cr is a KKT point 
of (P d ). If G(cr) y 0, then x = G^^Ficr) is the global minimizer of (V). 
If det(G(<x)) = 0, the global minimizer x of (P) is contained in the canonical 
equilibrium equation G(<r)x = F(&). 

Proof. By introducing Lagrange multiplier e e R" l+ " (where R m+ ™ is the 
nonpositive orthant of R m+ ") associated with cr > 0, the Lagrangian L(e,a) 
is given by 

L(e,cr) = -^F T (cr)G''(cr)F(cr)-cr T d-a-e T cr. (16) 

It is easy to prove that the criticality conditions V CT L(e, cr) — lead to 

ix T BiX - x - 6i \ 
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and the accompanying KKT conditions include 

< CT fe ± ^i T B k x - b£x - b k < 0, A; = 1, • • • , m + n. (18) 

From the complementary slackness, we can see that the x satisfies the con- 
straints in (V). Furthermore, since <r > for any g(x) < 0, we have 

P(x) > P(x) + <r T g(x) 

^ m-\-n _^ 

= 2 xTj4x - &Tyi ~ a + (2 xT * fe ' BfeX _ *fe b fc x _ &kh) 
fc=i 

= ix T G(a-)x - x t F(ct) - & T d - a 

= S( X ,a). (19) 

Noting that P(x) = S(x, <x), V x H(x, er) — and S'(x, er) is a quadratic func- 
tion with respect to x, we have 

P(x)-P(x) > «(x,cr)- »(x,<t) 

= (x - x)V x H(x, cr) + -(x - x) T V xx S"(x, <r)(x - x) 

= i(x-x) T G(<r)(x-x). (20) 

If G(<x) ^ 0, it is easy to find that x is the global minimizer of (P), where, x 
is contained in the canonical equilibrium equation 

G(ct)x = F(a). (21) 

If G(<x) is nonsingular, we have x = G~ 1 (cr)F(&). 



3 Implementation techniques 

To solve the optimization problem of (V d ), we firstly rewrite it to the following 
form: 

• 1 t , 
mm — t + <T d 

subject to : t > F T {ct)G- 1 {(t)F{ct) (22) 
G(cr) >r (23) 
er > (24) 

The global solution of (P ) is the same to the problem. Using the Schur 
complement [1], we can get the equivalent positive(semi) definite condition to 
([2"2]) and $T5j) 



G(<x) F(<t) 
F T (<t) t 



h (25) 
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and then the optimization problem can be expressed as the standard SDP 
form 

■ 1 T , 

mm —t + crd 

subject to : [St^ h (26) 

cr > (27) 

If G(a-) >~ 0, we can get the corresponding global solution to (V) by the 
canonical duality theory. In practice, the estimation of G(<r) may exist little 
inaccuracy due to the perturbed complementary slackness in primal-dual inte- 
rior point method and numerical precision. In this study, we use the Cholcsky 
factorization, Condition number and the smallest Eigenvalue of G(&) to eval- 
uate the positive defmitcness comprehensively. If G(<r) is ill conditioned or 
det(G(<x)) = 0, we can add a linear perturbation to the primal objective func- 
tion and then integrate the canonical dual solutions with the KKT conditions 
to recover the approximate solution or global solution to primal problem. De- 
tails of the techniques are given in the following examples. 



4 Numerical results 

All of the benchmark constrained optimization are from [9] , and we keep the 
number of each problem. In the experiments, we use SeDuMi [llj (a software 
package which can solve SDP problem) to obtain the canonical dual solutions. 
The built-in functions fsolve and fminunc in MATLAB are also used to solve 
the simple nonlinear equations and unconstrained optimization problems. 

Example 1 : gOl 

4 4 13 

min/(x) = x t - 5^x, 2 - ^x t . 

i—l i—1 z—5 



91 (x) 


= 2xi + 


2x 2 4 


- Xio 


f Xn 


- 10 < 


.92 (x) 


= 2xi + 


2x 3 4 


- Xio 


f X12 


- 10 < 


53 (x) 


= 2x 2 + 


2x 3 4 


- Xn 


f X12 


- 10 < 


54 (x) 


= -8xi 


+ Xio 


< o 






55 (x) 


= -8x 2 


+ %u 


< o 






56 (x) 


= -8x 3 


+ Xl2 


< o 






57 (x) 


= —2X4 


- x 5 - 


f Xio 


< o 




5s (x) 


= ~2x 6 


- x 7 - 


f- Xn 


< o 




59 (x) 


= -2x 8 


- x 9 - 


f Xl2 


< o 





where the bounds are < x 4 < l(i = 1, • • • , 9), < x. t < 100(i = 10, 11, 12) 
and < xi 3 < 1. 
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Solving the canonical dual problem, we can obtain <x = 



( al 


CT2 


0"3 


(7 4 


<?5 


^6 


<7 7 


08 


09 


(710 


(711 \ 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


1.0000 


1.0000 


1.0000 


5.0000 


5.0000 


012 


<7l3 


(Tl4 


015 


(716 


an 


<7l8 


019 


(720 


0"21 


(722 


\ 5.0000 


7.0001 


2.0001 


3.0001 


2.0001 


3.0001 


2.0001 


-0.0000 


-0.0000 


-0.0000 


1.0001 / 



In this case, G(&) >z but singular, satisfying the global optimality condi- 
tion. By the KKT condition, we can find that 37,58,59, bounds of xi,- ■ ■ ,xg, 
and X13 are active, so we can first get 



( Xi 


Xl 


%3 


Xa 


x 5 


x 6 


x 7 


#8 


Xq 


£io 




«12 


«13 \ 




1 


1 


1 


1 


1 


1 


1 


1 








1 J 



where, " — " means undetermined. Considering that constraints 37,38,39 are 
active, solving the corresponding linear equations, we can easily get Xio = 
3,Xn — 3,x\2 — 3. Finally, the global solution to gOl is x* = 



( X\ 


Xi 


X3 


X4 


x 5 


X6 


X7 


Xg 


£9 


£10 


Xu 


£12 


Xl3 \ 




1 


1 


1 


1 


1 


1 


1 


1 


3 


3 


3 


1 y 



and /(x*) = -15. 

Example 2: g04 
min/(x) = 5.3578547x^ 



subject to : 31 (x 

32 (x 

33 (x 

34 (x 

35 (x 

36 (x 



0.8356891x1X5 + 37.293239xi - 40792.141 
85.334407+ 0.0056858x 2 x 5 + 0.0006262xix 4 - 0.0022053x 3 x 5 - 92 < 
-85.334407- 0.0056858x 2 x 5 - 0.0006262xix 4 + 0.0022053x 3 x 5 < 
80.51249 + 0.0071317x 2 x 5 + 0.0029955xix 2 + 0.0021813X 2 . - 110 < 
-80.51249 - 0.0071317x 2 x 5 - 0.0029955xix 2 - 0.0021813x2 + 90 < 
9.30096 + 0.0047026x 3 x 5 + 0.0012547xix 3 + 0.0019085x 3 x 4 - 25 < 
-9.30096 - 0.0047026x 3 x 5 - 0.0012547xix 3 - 0.0019085x 3 x 4 + 20 < 



where 78 < x x < 102, 33 < x 2 < 45 and 27 < x l < 45(z = 3, 4, 5). 
Solving the canonical dual problem, we can obtain & = 



( G1 


(72 


0"3 


(7 4 


05 


0"6 


(7 7 


0"8 


09 


cio 


(711 \ 


^336.8388 


0.0000 


0.0001 


0.0002 


0.0003 


798.2826 


2.0310 


6.1233 


0.0001 


1.6054 


1.1849 J 



In this case, G(er) >- and cond(G(er)) = 9.7330e5, satisfying the global 
optimality condition, so we can get x = 



( Xl 


Xl 


X3 


X4 


x 5 \ 


^77.9452 


33.0179 


29.7345 


44.9884 


38.2523 J 



Noting that the condition number is large, according to the KKT condition, 
we can first get 



( Xl 


Xl 


X3 


Xi 




\78 


33 




45 


X5^ 
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Considering that constraints g\ , g§ are active, solving the corresponding linear 
equations, we can easily get x 3 = 29.995256025681599, x 5 = 36.775812905788207. 
Finally, the global solution to g04 is x* = 



( Xl 


Xi 


\78 


33 



and /(x*) 



X3 

29.995256025681599 
3.0666e4. 



x 5 

36.775812905788207 



Remark 1 We use the inverse of G(er) because only its smallest eigenvalues 
approximates to zero although its condition number is large. As a matter of 
fact, the solution x causes only little infeasibility of the first constraint. By 
integrating the canonical dual solutions and the KKT conditions, we claim 
that Xi,X2 and X4 are determined in the first stage. 

Example 3: g07 

min/(x) = xj + x\ + x x x 2 - 14xi - 16x3 + (x 3 - 10) 2 + 4(x 4 - 5) 2 + (x 5 - 3) 2 
+ 2(x 6 - l) 2 + 5x? + 7(x 8 - ll) 2 + 2(x 9 - 10) 2 + (bio - 7) 2 + 45 
= -105 + 4xi + 5x 2 - 3x 7 + 9x 8 < 
= 10xi - 8x 2 - 17x 7 + 2x 8 < 
= -8x1 + 2x 2 + 5x 9 - 2xio - 12 < 
= 3(xi - 2) 2 + 4(x 2 - 3) 2 + 2xj - 7x 4 - 120 < 
= 5x 2 + 8x 2 + (x 3 - 6) 2 - 2x 4 - 40 < 
= x 2 + 2(x 2 - 2) 2 - 2xix 2 + 14x 5 - 6x 6 < 
= 0.5(xi - 8) 2 + 2(x 2 - 4) 2 + 3X 2 , - x 6 - 30 < 

2 - 7xm < 



subject to : g\ (x 

92 (x 

93 (x 

94 (x 

95 (x 
9e (x 
97 (x 
9s (x 



= — 3xi + 6x 2 + 12(xg 



where -10 < x; < 10(i = 1, • • • , 10). 

Solving the canonical dual problem, we can obtain <x = 



/ 01 


ct 2 


^3 


<T 4 


CT5 


^6 


(7 7 


08 


0-9 


\ 


1.7168 


0.4746 


1.3760 


0.0205 


0.3120 


0.2871 


0.0000 


0.0000 


0.0000 




O"10 


CTll 


C12 


Cl3 


(7i4 


CT15 


0"16 


(717 


C18 




\ 0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


/ 



In this case, G(<x) >- and cond(G(er)) = 7.0000, satisfying the global 
optimality condition, so we can get x* = 



( Xi 


X2 


X3 


X4 


x 5 


x 6 


x 7 


x 8 


x 9 


XlO \ 


^2.1721 


2.3636 


8.7746 


5.0959 


0.9903 


1.4307 


1.3218 


9.8286 


8.2800 


8.3760 J 



and /(x*) = 24.3111. Note that there exists little infeasibility due to numerical 
precision. 
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Example 4 '■ glO 



min /(x) = x\ + x 2 + x 3 
subject to : gi(x) = -1 + 0.0025(x 4 + x 6 ) < 

g 2 (x) = -1 + 0.0025(x 5 + x 7 - x 4 ) < 
3 (x) = -l + 0.01(x 8 -x 5 ) <0 

,g 4 (x) = -X!X 6 + 833.33252x 4 + lOOxi - 83333.333 < 
g 5 (x) = —x 2 X7 + 1250x5 + x 2 x 4 — 1250x 4 < 
g 6 (x) = -x 3 x$ + 1250000 + x 3 x 5 - 2500x 5 < 

where 100 < x x < 10000, 1000 < x, < 10000(i = 2, 3) and 10 < x t < 1000(i = 
4,--- ,8) 

Solving the canonical dual problem, we can obtain a = 



( *i 
9.2834 

08 

\0.0001 



02 
28.9205 

0"9 
0.0001 



0"3 

5.8893 

01o 
0.0000 



(74 

0.0001 

o n 
0.0000 



05 

0.0001 

012 

0.0000 



0"6 

0.0001 

013 

0.0000 



cry 

0.0001 

014 

0.0000[/ 



In this case, G(&) >- and cond(G(<r)) = 749.4514, satisfying the global 
optimality condition. However, the max(eig(G(<x))) = 2.5743e — 4, which is too 
small, so we cannot use the inverse of G{&) directly. By the KKT condition, 
we can find that constraints gi, g 2l g 3 are active, and all of the box constraints 
are inactive. That is to say, the problem is equivalent to a linear programming 
problem with linear constraints, which indicates that #4, <75, ge must be active. 
Fixing X4 , X5 , we have 



Xi = 

x 2 = 

x 3 = 
x 6 = 
x 7 = 
x 8 = 



83333.333 - 833.33252x 4 

x 4 - 300 
1250x 4 - 1250x 5 

x 5 - 400 
12500 - 25x 5 
400 - x 4 
400 + x 4 - x 5 
100 + x 5 



As a result, we can reduce the problem to 

. ,, , 83333.333- 833.33252x4 1250x4 - 1250x5 
mm/(x) = 1 — h 12500- 25x 5 



x 4 - 300 



x 5 - 400 



Taking the box constraints of x\ , • • • , x% into consideration, when using (100, 200) 
as an initial point for the unconstrained optimization problem with two vari- 
ables, it is easy to get the only minimum x 4 = 182.0176995811199 and x 5 = 
295.6011732779338. Utilizing the equations obtained by the complementary 
slackness, finally, we have x x = 579.3066844253549 x 2 = 1359.970668051655, 
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x 3 = 5109.970668051655, x 6 = 217.9823004188801, x 7 = 286.4165263031861, 
x s = 395.6011732779338 and /(x*) = 7049.248020528666. 

Remark 2 Wc don't use the inverse of G(<x) directly because all of its eigen- 
values are approximately zeros. And the reason why we still use the canonical 
dual solutions as useful heuristics is that the G(&) is slightly positive definite 
due to the perturbed complementary slackness caused by the SeDuMi. Since 
all of the box constraints are inactive and the objective function is linear, it 
is not difficult to imagine that all of the constraints must be active. Note that 
the constraints of a: 4 and x^ are changed when solving the unconstrained opti- 
mization problem since constraints of x\, X2, x 3 and xq, x 7 , x% must be satisfied. 

Example 5: gl8 

min/(x) = — 0.5(a;ia;4 — x 2 x 3 + x 3 x 9 — x 5 x 9 + x 5 x$ — x e x 7 ) 
subject to : 31 (x) = x\ + x\ — 1 < 
52 (x) = x\ - 1 < 
g 3 (x)=x 2 5 +x 2 6 -l<0 
,g 4 (x) = x\ + (x 2 - x 9 ) 2 - 1 < 
3 5 (x) = (xi - x 5 ) 2 + (x 2 - x 6 f - 1 < 
36 W = {xi - x 7 ) 2 + (x 2 - x$) 2 - 1 < 
g 7 (x) = (x 3 - x 5 ) 2 + (x 4 - x 6 ) 2 - 1 < 
g 8 (x) = (x 3 - x 7 ) 2 + (x 4 - x$) 2 - 1 < 
3 9 (x) = x 2 + (x s - x 9 ) 2 - 1 < 
.gio (x) = x 2 x 3 - xix t < 
3n (x) = -x 3 x 9 < 

312 (x) = x 5 x 9 < 

313 (x) = x 6 x 7 - x 5 x$ < 



where -10 < x-i < 10, (i = 1, • • • , 8) and < x 9 < 20. 

Solving the canonical dual problem, we can obtain & = 



( al 


CT2 


^3 


C7 4 


^5 


^6 


cr 7 


^8 


09 


010 


011 \ 


0.1444 


0.0000 


0.1444 


0.1445 


0.0000 


0.1442 


0.1441 


0.0000 


0.1445 


0.0000 


0.0000 


012 


<7l3 


(Ti4 


Cl5 


Cl6 


an 


<?18 


019 


020 


021 


022 


\ 0.0000 


0.0000 


-0.0000 


-0.0000 


-0.0000 


-0.0000 


-0.0000 


-0.0000 


-0.0000 


-0.0000 


0.0000 / 



In this case, G(&) >- and cond(G(er)) = 7.1887e7, satisfying the global 
optimality condition. However, the condition number is large. Taking the KKT 
conditions into account, we can conclude that constraints 31, 33, 34, ga, g 7 , 
g 9 are active since the corresponding a\, <j 3 , 04, <j 7 , a 9 are not zeros. But 
it becomes still difficult to solve the nonlinear equations. Considering that 
several eigenvalues of G{&) are zeros and there exists no linear term in the 
objective function, and in this situation, we add a small linear perturbation 
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0.05 (xi + ••■ ,xq) to the primal objective function. Solving the perturbed 
canonical dual problem, we get G(<x) >- and cond(G(er)) = 1.4592e3 and 
the smallest eigenvalue of G(<x) is 0.0021. Therefore, we can get x = 



f Xl 


X 2 


%3 


X4 


x 5 


x 6 


x 7 


x$ 


Xg 




^-0.9660 


-0.2585 


-0.2587 


-0.9660 


-0.9661 


-0.2588 


-0.2589 


-0.9657 


0.0005 





and /(x) = —0.8663. Note that there exists little infeasibility due to numerical 
precision. 



Remark 3 The solution we get is quite different from the best known solution. 
According to the canonical duality theory, the global solution to this problem 
is not unique. The linear perturbation technique can only help to find one of 
the global solutions. As a matter of fact, the following solutions x = 



( Xx 


X 2 


xz 


X 4 


x b 


x 6 


x 7 


x s 


Xg 




\ 0.0450 


-0.0387 


0.8663 


-0.4999 


0.0004 


-1.0001 


0.8878 


0.5000 


0.9604 





( Xl 


X2 


X3 


X4 


x§ 


x 6 


x 7 


x & 


Xg 




^0.0689 


-0.9972 


0.9088 


-0.4179 


0.0920 


-0.9959 


0.8986 


-0.4388 


0.0009 


) 



and 



( Xl 


X2 


X3 


Xi 


x 5 


x 6 


x 7 


x 8 


Xg 




^0.6888 


-0.7257 


0.9693 


0.2454 


0.6973 


-0.7173 


0.9726 


0.2332 


-0.0006 


) 



can all be considered as approximate solutions, which are obtained by the 
proposed techniques. 

5 Conclusion 

We have applied the canonical duality theory to solve a class of CEC bench- 
mark constrained optimization problems. Experimental results show that some 
of the examples can be solved directly, some of them can be solved by inte- 
grating the canonical dual solutions and KKT conditions, and some can be 
solved approximately by adding a small linear perturbation term. 
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